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Elastic constants and zone-boundary phonon frequencies of gold are calculated by total energy electronic 
structure methods to twofold compression. A generalized force constant model is used to interpolate throughout 
the Brillouin zone and evaluate moments of the phonon distribution. The moments are used to calculate the 
volume dependence of the Griineisen parameter in the fee solid. Using these results with ultrasonic and shock 
data, we formulate the complete free energy for solid Au. This free energy is given as a set of closed form 
expressions, which are valid to compressions of at least V/Vb = 0.65 and temperatures up to melting. Beyond 
this density, the Hugoniot enters the solid-liquid mixed phase region. Effects of shock melting on the Hugoniot 
are discussed within an approximate model. We compare with proposed standards for the equation of state to 
pressures of ~ 200 GPa. Our result for the room temperature isotherm is in very good agreement with an earlier 
standard of Heinz and Jeanloz. 

I. INTRODUCTION 

The elastic constants, phonon frequencies, and equation of state (EOS) are fundamental properties of matter. The values 
of these parameters under compression find application in geophysics 1 and in the prediction and interpretation of processes 
of dynamic compression. 2,3 For many materials, especially elemental metals, the principal Hugoniot curves, the set of states 
accessible via a single shock from ambient conditions, have been measuredi 2 ^^ Since the pressure, density, and internal energy 
are known along the Hugoniot from the jump conditions, 2 this data is an important baseline for high pressure equations of state. 
The off-Hugoniot EOS is needed for the prediction of processes involving more complicated loading paths, such as multiple 
shocks or shock and release. 7 Extrapolation from the Hugoniot has also been used to establish pressure standards for static high 
pressure experiments, whose room temperature isotherms are regarded as known. 2 - 8 - 9 Au has been used as a standard in this 
way, and has been used to calibrate the ruby Ri line as a secondary standard. 91011 Recent studies have called into question the 
accuracy of the gold pressure standard. Akahama et alM> compressed Au and Pt in the same cell and found two Au standards!^ 
to give pressures lower than Pt by 20 and 15 GPA, respectively, at 150 GPa. Shim et aLp^on the other hand, propose a new EOS 
for Au that gives pressures still lower than either of these earlier standards. 

To relate the Hugoniot to the room temperature isotherm requires information on the Griineisen parameter y = V (dP/dE) v 
and the specific heat Cy? These are dominated by lattice vibrations in the regimes under consideration here. With increasing 
compression, the separation between the Hugoniot pressure and the room temperature pressure increases, resulting in greater de- 
pendence of the inferred room temperature isotherm on y. Because y is not easily measured at high pressure, this introduces a non- 
negligible source of uncertainty in the pressure standards. In most cases, a model assumption of the form y(V) = y(Vo)(V/Vo) l? 
has been used. The specific value q = 1 has been used often in shock worki 2 ^ However, a power law dependence of y on V is 
qualitatively incorrect at high compression, and extrapolation on this basis is inherently limited. 

In principle, the various components of the EOS - the static lattice energy, the vibrational free energy, and the electronic 
excitation free energy - can be evaluated from electronic structure theory. Practical calculations based on approximate density 
functional theories typically have errors of several percent in the density at zero pressure, and 10 % or more in the bulk modulus. 15 
These are unacceptably large errors for the purpose of high accuracy equations of state. These properties can be measured 
accurately, and are mainly determined by the cold energy curve, which is the largest contribution to the EOS in the regime 
considered here. On the other hand, we find that ab initio electronic structure calculations are capable of obtaining phonon 
frequencies of sufficient accuracy to strongly constrain the volume dependence of the Griineisen parameter y. We therefore 
propose that the most accurate EOS in the solid is obtained by combining an empirical cold energy with ab initio lattice vibration 
and electronic excitation free energies. This is analogous to procedures that have been used to reduce shock Hugoniot data and 
derive room temperature standards, but the present analysis gives a strong physical foundation for the volume dependence of y, 
allowing confidence in the results at higher compression. We are not aware of any evidence of solid-solid phase transitions in 
Au, and we consider only the fee solid. Recent calculations 16 indicate a transition to the hep structure for V/Vq < 0.6, and we 
do not attempt to extend our solid EOS beyond this density. 

This paper is arranged as follows. First we discuss lattice dynamics and its connection to the EOS. We emphasize the 
importance of the classical limit for defining the ion motion contribution to y. We describe our procedures for calculating 
phonon frequencies and interpolating to the whole Brillouin zone. Next we present our results for elastic and phonon properties 
of Au, and our analysis to obtain y lon (y). We then present the complete EOS for solid Au by giving a set of closed form 
expressions for the Helmholtz free energy with parameters. The resulting room temperature isotherm is compared with various 
proposed standards. Next we discuss shock melting and its effect on the Hugoniot within an approximate model, and finally give 
our conclusions. 
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II. LATTICE DYNAMICS AND THE EQUATION OF STATE 



We write the Helmholtz free energy as 



F(V, T) = (j)o(V) + F ion (VJ) + F el (V,T) 



(1) 



where cbo is the static lattice energy, F lon is the ion motion free energy, and F el is the electronic excitation free energy. In 
our applications, (|)o gives the largest contribution to the pressure. The ion motion free energy gives the dominant temperature 
dependence of P — — (dF/dV) T . The electronic excitation term F ei is generally a small correction to the solid EOS, becoming 
important for the liquid Hugoniot at several hundred GPa. 

It is, therefore, imperative to have an accurate F 1011 for calculating the temperature dependence of P. The quasi-harmonic 
approximation has been found to have small errors in many cases. 17 We have carried out Monte Carlo simulations using an 
embedded atom model 18 of Cu to investigate the importance of anharmonic corrections to the ion free energy. Here we define 
the anharmonic free energy to be the difference between the true ion free energy and the quasi-harmonic approximation, noting 
that the term has been used differently 19 by other authors. 20 We intend to publish details of our Cu simulations elsewhere. In 
summary, we find that up to 300 GPa on the melting curve, |f anh | < 0.85 GPa, and is never more than 3% of the thermal pressure 
P(V, T) — P(V, 0), a very small correction to the total pressure. Given the similarity of the bonding in Cu and Au, we expect these 
results to be relevant for Au also. Thus, in what follows we neglect anharmonicity and use the quasi-harmonic approximation 
forF 10 ", 
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where we have introduced g((o) — j^Y^k S(<*> — K>k) , the normalized phonon density of states. The phonon frequencies cot are 
functions of volume only, and the sum is over the 3N normal modes of the crystal. 

Application of Eq. (0 requires the full phonon density of states g((o) at all volumes. While in principle this information is 
available from our calculations, in practice the classical limit dominates our EOS, allowing for a substantial simplification. In 
the classical limit, which is the leading term in the high temperature expansion of Eq. 0, the free energy is given by 



F cl (V,T) = 3Nk B Tln 



- l/3 W 

k B T 



(3) 



where we have introduced the moment 
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Other moments are conventionally defined as 
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(5) 



where the normalizations in Eqs. l|4} and (0 are chosen so that for a Debye spectrum, described by 



g((0) = -rC0 2 ©(C0 D -C0), 



(6) 



all co„ are equal to C0d, the Debye frequency. We subsequently give results for v„ = C0„/ (2k) corresponding to the experimental 
convention of giving frequencies v in Hertz, as opposed to angular frequencies, CO. 

In the classical limit, the ion pressure is linear in T with (dP ion /dT) v = (3Nk B /V)y , where y = ~^ L , and the specific 
heat is constant, Cy = 3Nk B . The role of quantum ion motion in the EOS is illustrated in Figure^ which shows the pressure 
along an isochore V/Vb = 0.8. The solid curve is our full EOS, described below, and the dashed curve uses the classical limit 
of F lon at all T. The arrows mark room temperature and the Hugoniot temperature, which is 1340 K at the given density. The 
melting temperature at this density is estimated to be 4900 K, beyond the range of the plot. It is clear that the classical limit 
dominates even at room temperature. The Hugoniot is well into the classical regime. At higher densities, the Hugoniot is still 
further above the quantum regime. The largest quantum effect on the pressure is at T = 0, where the zero-point vibrations 
contribute a pressure of 0.8 GPa at this density. The classical limit is clearly dominant for interpolating between the Hugoniot 
and room temperature. 
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FIG. 1: Role of quantum ion motion in the EOS. Graph shows P(T) along an isochore V /Vq = 0.8. Solid curve is the full EOS, dashed curve 
uses the classical limit of the ion free energy. Hugoniot temperature at the given density is 1340 K. 



We have found that an interpolation based on the Debye model gives a very accurate approximation to the full quasi-harmonic 
free energy. The Debye free energy is a special case of the quasi-harmonic free energy, Eq. (0. Inserting the Debye density of 
states, Eq. (|6}, gives 



F D (V,T) =N 

where 



h(£> D + 3k B T\n [ 1 - e -'" s>DlkBT ) - k B TD (h& D /k B T) 



(7) 
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In light of the above remarks, it is important to capture the correct classical limit. This requires that we set 

co D (V) = co (V). CO 

To emphasize the distinction between Eq. © and the standard definition of cod in terms of acoustic velocities, we refer to the 
Debye free energy together with Eq. (|9) as the high-temperature Debye model. The high-temperature Debye model gives the 
same results as the full quasi-harmonic free energy in the classical regime, and in addition obeys the Nernst theorem at low T . 
We have checked that at low T, the error in the pressure compared to the full quasi-harmonic approximation is entirely negligible. 

Thus, we can simplify the specification of the lattice vibration free energy from giving g(co) at all V to giving the single 
moment coo at all V. This is advantageous for numerical applications. It also allows us to express our full EOS in a few compact 
formulae so that it is generally accessible. These formulae are given below. 

Calculating the moments co„ requires the phonon frequencies for all k in the Brillouin zone. Direct ab initio calculations on 
a dense mesh in k would be quite expensive. As a result of another study 23 we found that the low order moments can often be 
accurately obtained from short-ranged force constant models. In particular, for Au, the moment ©o is converged to less than 1 % 
with a 2nd neighbor (2NN) model. Therefore, our method is to calculate four zone boundary phonon frequencies corresponding 
to the transverse and longitudinal modes at the X and L points. These are calculated with standard frozen-phonon methods. In 
addition the three elastic moduli are calculated using the method described by Soderlind et al>2i We fit these results to a 2NN 
force constant model, which then allows the evaluation of <Mu for arbitrary k for integration over the Brillouin zone. 

The electronic structure calculations used the full-potential LAPW code WIEN97. 25 We used the LDA rather than the GGA, 
based on Boettger's finding 16 that the LDA gives a better static lattice energy than the GGA for Au. Some numerical parameters 
used in the calculations were, in atomic units: muffin tin radius tmt = 2.0; plane wave cut-off oyiT^max = 9.0; cut-off for 
expansion of density and potential g max = 14. For elastic modulus calculations, Brillouin zone integrals used special points 
corresponding to 18 3 points in the full zone, with Gaussian smearing of the energies by 20 mRy. The zone boundary phonons 
were found to be comparatively insensitive to the k-point mesh, and smaller meshes of 10 3 points were used. The 5/?, 5c?, and 
6s shells were treated as valence states, and local orbital extensions 26 were used in the p and d channels. 



4 




FIG. 2: (Color Online) Phonon dispersion for Au at room temperature and ambient pressure. Open circles are experimental data of Lynn et. 
almL Filled diamonds are LDA calculations of zone boundary phonons at the X and L points. Solid curve is interpolation based on fit of 2NN 
force constant model to the LDA zone-boundary phonons and elastic moduli. 



We used a generalized Born-von Karman force model for our lattice-dynamical calculations, from which we were able to 
compute the phonon dispersions in the entire Brillouin zone. Since gold has a very simple phonon dispersion, we employed only 
a model with first (INN) and second (2NN) nearest-neighbor interatomic shells of atoms. In an fee lattice the INN and 2NN 
forces are determined by three and two independent parameters, respectively (for more details see e.g. ref.l23l). These five force 
constants were extracted by fitting simultaneously the phonon frequencies of gold at the X and L points of the Brillouin zone and 
the elastic constants Cn, C12, and C44 near the zone center (r point). Thus we fit a total of seven independent data points. We 
gave equal weight to the zone boundary phonons and elastic constants in our % 2 -fit of the phonon dispersions. Although the 2NN 
Born-von Karman force model is too simple to reproduce all phonon frequencies within less than approximately 10%, see Figure 
12 it is sufficiently accurate to compute integrated quantities such as the phonon moments within less than approximately 3%. 
More accurate phonon dispersions can be obtained if needed, by computing ab initio frequencies at half and quarter distances 
in the Brillouin zone and fitting those to a 3NN or higher-order Born-von Karman force model, or to phonon models with 
interatomic pseudo-potentials. 



III. ELASTIC AND PHONON PROPERTIES 

Figure |2 shows the calculated phonon dispersion curves for Au at the density corresponding to ambient pressure and room 
temperature. The filled diamonds at the X and L points are the LDA frozen phonon results. These, together with the calculated 
elastic moduli are used to obtain the force constants. The solid curves are the interpolation to general wavevectors using the 
force constant model. The open circles are the experimental data of Lynn et a/j2lThe force constant model is over-constrained, 
so it neither goes precisely through the calculated zone boundary frequencies, nor has it exactly the calculated elastic moduli. 
The shapes of the dispersion curves are simple enough for Au that they are generally well captured by the 2NN model. 

Table U summarizes our results for the zone boundary phonons, elastic moduli, and moments v„ = co„/2ti, as functions 
of volume. The reference volume Vo corresponds to T = 298 K and P = 1 bar = 100 kPa. For Au, Vo = 10.212 cm 3 /mol 
or 1 14.43 flg/atom. Also shown are the experimental data at room temperature. No attempt is made to account for the 
temperature dependence beyond comparing at the correct volume. There is generally good agreement between the calculated 
and experimental quantities. The main exception is C44, which is calculated substantially lower than measured. Our calculation 
is in better agreement with an earlier LDA calculation. 24 The main result of the present calculation is the value of Vo, which 
is within 3% of the measurement. The ratio Vi/vo is 1 for a Debye spectrum. Our calculations give Vj/Vo = 1.03 at ambient 
density and 1 .04 at twofold compression, so, by this measure, the departure from a Debye spectrum is small and nearly constant 
with volume. 

Measurements of elastic moduli were reported by Duffy et al.M> to P = 37 GPa, and ab initio calculations by Tsuchiya 
and Kawamura 30 to 100 GPa. Figure [3] shows these results along with our calculations. Our volumes are converted to room 
temperature pressures using the EOS described below. The highest density in table [Q corresponds to P « 780 GPa, and the 
graph is restricted to lower pressures to highlight the comparison with these other works. Our calculations of C\\ and C12 
are in good agreement with Tsuchiya and Kawamura, while our C44 is systematically lower, in somewhat better agreement 
with the experiments. Our C\\ is in good agreement with the experiments while C12 is slightly high. It should be noted that the 
experimental C, ; - depend on a model parameter a used in the analysis. 29 The plotted points correspond to a — 1 . It is encouraging 
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TABLE I: Calculated phonon and elastic properties of Au. Vo is the volume at 298 K and atmospheric pressure and is 10.212 cm /mol for Au. 
Frequencies v given in THz and elastic moduli given in GPa. 
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FIG. 3: Pressure dependence of the elastic moduli of Au. Open squares are present LDA calculations with pressures from our EOS. Solid 
circles are experimental data from Duffy et al*^. Open triangles are previous LDA calculations by Tsuchiya et al. 3() 



that our calculations give the correct trends for the pressure dependence of Cy. 

Having determined the moments v„ as functions of volume, we need to interpolate and differentiate them to obtain thermody- 
namic functions. To do this, we assumed the following functional form for J„(V), which has been used in creating many wide 
ranging equations of stated 



Jn(V)=f 



B„ 



(10) 



where "f is the infinite density limit of y, and A„ and B„ are parameters. Alternatively, we can express A„ and B„ in terms of 
q„ = <iln(y„)/flfln(V) andy,(V ) as 



A n = y„(Vo)[2-<?„(yo)]-2y 00 

B„ = y n (Vo)[qn(Vo)-l] + T 



Integrating y„ 



d In V, 
" c/lnV 



f, we have 



v„(V)=v„(V ) 



-r 



exp < 



2 



(11) 



(12) 



By fitting this functional form to the calculated Vo(V), we extract the parameters giving yo(V ). In our fits, we keep y°° fixed. The 
value y 00 — 2/3 has been widely used, 32 although arguments have been made for y°° = 1/2. 33 For our application and density 
range, we have found that the quality of the fit and the resulting yo(V) are insensitive to this choice, and we use y°° = 2/3. The 
fitting procedure is illustrated in Figure |4j where we show the fit for Vo(V), the log moment and the resulting yo(V). We did 
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FIG. 4: Volume dependence of Vo, the log moment of the phonon frequencies, and Yo = — rflnVo/rflnV for Au. Lower graph is Vn(V). Open 
symbols are calculated moments. Upper graph is resulting Yo(V). Solid curves are present results fitting Eq. I lOt to the calculated moments 
with the constraint Yo(Vo) = 2.97. Dash-dot curve shows the commonly used approximation j(V) = j(Vo)(V /Va) q for q = 1.7. 13 Dashed curve 
corresponds to the y of Jamieson et al&, who used q = 1 and j(Vq) =3.215. 



the fits with Vo(Vo), Yo(Vo), and qo(Vo) as free parameters, and with Yo(Vb) constrained to the experimental value of 2.97. (This 
experimental value was determined so that the thermal expansion for the subsequent EOS overlies the recommended curve of 
Touloukiar. 34 from 100 to 1200 K.) The unconstrained fit gives Yo(Vb) = 2.88. Constraining Yo(Vo) increases the RMS error of 
the fit from 4.9 x 10~ 2 to 5.7 x 10~ 2 THz, which is not a significant increase. Also, we note that constraining Yo(Vo) does not 
change Yo(V) at smaller V. For this reason, we believe that the constrained fit is the best overall approximation for Y° n (V) up 
to twofold compression. This is shown as the solid curve in Figure @] The corresponding parameters are Vo(Vo) = 3.46 THz, 
Jo(Vo) = 2.97, go(Vo) = 1 .37. We do not show the unconstrained fit in Figure|4]because it is too close to the constrained fit. The 
dot-dashed curves correspond to q = const, y(V) = j(Vo)(V /Vo) q . The value q = 1 gives a y that is too high over this density 
range, while q — 1.7 13 is too low. Jamieson et alfi used q = 1 and an ambient value y(Vb) = 3.215, which results in the dashed 
curve. We believe that their y is too high at all volumes considered here. Any constant value of q will lead to a y that is too small 
at very high compression. 



IV. EQUATION OF STATE 

The complete EOS is determined when the Helmholtz free energy F is given as a function of volume and temperature. In this 
section we describe our form for F and give the numerical parameters for fee Au. First we write the free energy as in Eq. Q as 
the sum of the static lattice energy (])o, an ion motion free energy F lon , and an electronic excitation term F el . For the ion motion 
free energy we use the quasi-harmonic approximation, and further specialize to the high temperature Debye model. Thus we 
take the ion free energy to be given by Eqs. Q and (|8j, where we identify C0d with the log moment ©o of the phonon frequencies. 
As discussed in section II, this gives the correct classical limit, and leads to extremely small errors at low T compared to the 
full phonon spectrum. This way we are able to give a closed form expression for F lon , which is very convenient in numerical 
applications. The volume dependence of coo is that given in Eq. i ll 2t with parameters given above. For the final EOS we have 
replaced the fitted value of (Mo(Vo) with the experimental one. This is essentially a small shift in the absolute entropy, which has 
negligible effect on the results. 
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Electronic excitations give a small contribution to the thermodynamics, which we approximate by 

F e \v,T) = -^Nr(V)T 2 . (13) 

The Sommerfeld coefficient T is proportional to the electronic density of states at the Fermi energy, T = (% 2 /3)&|g(E/). Our 
calculations show that dlng(£f)/dlnV « 0.76 over the range of densities considered, so we take 

r(v) = r(v )(v/v ) K (14) 

with K = 0.76. The calculations give r(Vb) = 6.7 x 10~ 4 J/mol K 2 , while the measured value from the low-temperature specific 
heat is— r(Vb) = 7.28 x 10~ 4 J/mol K 2 . The measured low-temperature specific heat coefficient r expt is expected to be enhanced 
with respect to the bare value T as a result of electron-phonon interactions by a factor (1 + X), where X is the dimensionless 
electron-phonon mass enhancement parameter. The value for gold is X « 0.05 — 0.15, 36 which agrees with our calculated 
r/r expt = 1.09. At high temperatures, kgT > h&o, the phonon mass enhancement becomes ineffective. Hence we use our 
calculated density of states in the remainder of this paper. 

The static lattice energy is needed to complete the free energy. We adopt the functional form due to Vinet et al? 1 



X = \{B\-\) 



,y/3 

V" 



(15) 



which is parameterized by the volume at the minimum, V*, the bulk modulus, B* and its pressure derivative, B\. These pa- 
rameters have been determined empirically by requiring that the EOS reproduce the ambient volume, and ultrasonic data for 
B$ at ambient conditions. 28 For the latter, we adopt the value B$ = 173 GPa. Ultrasonic data gives values for (dB$/dP) T from 
5.2 to 6.4. Therefore, to complete the EOS, we require that it reproduce the measured slope of the Hugoniot, which we take 
from the fisi U s = 3.12km/s+ 1.521f/ p relating the shock velocity U s to the particle velocity U p . The resulting parameters give 
(dBs/dP) T = 5.49, consistent with ultrasonic data. 

The complete set of parameters for the free energy of fee Au are summarized here. The static lattice energy is given by 
Eq. (JBJ with 

V* = 10.0834 cm 3 /mol 
B* = 180.0 GPa 

B\ = 5.55. (16) 

The ion motion free energy is given by Eqs. and (jSJl with the volume dependence of C0d given by Eq. d!2l > (with Vo — ► COd, 
Yo — » y, etc.). The parameters A and B are given in terms of y(Vo) an d q(Vo) by Eq. il It . The numerical values are 

V = 10.212 cm 3 /mol 
co£,(Vo) = 22.9 x 10 12 s-' 
j(V ) = 2.97 
q{V ) - 1.3677 

Y°° - 2/3. (17) 
Finally, the electronic excitation free energy is given by Eqs. Jl 3i and dl4> with 

r(Vb) = 6.7 x 10" 4 J/mol K 2 

k = 0.76. (18) 



The following calculations use these parameters for the free energy of fee Au. Once the free energy is known, the pressure, 
internal energy, etc., can be evaluated. We calculate the Hugoniot by fixing a value of the volume and solving for the temperature 
such that the jump condition, E — Eq — \ (P +Pq) (Vo — V) > is solved. Here Eq, Vq and Pq correspond to the initial state, taken to 
be ambient temperature and pressure. 

Figure |3 shows the Hugoniot and room temperature isotherm for the present EOS. The solid symbols are the Hugoniot 
data. 5 - 6 ' 38 The dashed line corresponds to the linear fit for U S (U P ). Also shown as the open circles are the room temperature data 
from Heinz and Jeanloz 13 taken to 70 GPa with the ruby pressure scale. The open squares are the room temperature data of Bell 
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FIG. 5: Hugoniot and room temperature isotherm of Au. Solid curves are present solid EOS. Solid symbols are Hugoniot data^i— Note that 
there are five overlapping Hugoniot points at V /Vq « 0.7. Open symbols are room temperature dataiUJ— For the Bell et al, data, 11 pressures 
are from the extrapolated ruby standard. The onset of shock melting is estimated at P ~ 280 GPa, as marked. 



et al. n to 188 GPa. The pressure shown in the graph is taken from the ruby scale, which is extrapolated from a lower pressure 
calibration. 9 Our calculations agree with these data to 3% at the highest pressure, indicating support for the extrapolation of the 
ruby scale. 

For V > 0.656Vo, the Hugoniot lies in the solid. By matching the measured Hugoniot with our theoretically based Griineisen 
parameter, we confirm the accuracy of our room temperature isotherm, and we have high confidence in our EOS to this density. 
Combining our calculated Hugoniot temperatures with a Lindemann melting curve, we estimate that the Hugoniot enters a solid- 
liquid coexistence region at V /Vq = 0.656, P = 280 GPa. A more extensive discussion of shock melting is given in the following 
section. There we find that explicitly accounting for melting leads to good agreement with the high pressure Hugoniot data, 
indicating that our solid EOS is valid to densities of V /Vq « 0.6. 

Figure [6] shows our room temperature isotherm along with data from Heinz and Jeanloz, 13 and Bell et al*^ Also shown are 
some of the proposed EOS standards. Jamieson et al. & used the Hugoniot as a reference, and calculated the room temperature 
isotherm using a Mie-Griineisen EOS. They assumed j/V = const, and used y(Vb) = 3.215, so their y is always larger than 
ours. This results in lower room temperature pressures than ours, but by restricting their analysis to V/Vb > 0.775, the impact 
of the assumed y is minimized. The Heinz and Jeanloz EOS is based on extrapolating their room temperature data, with some 
consideration of the shock data. Our room temperature isotherm is in good agreement with Heinz and Jeanloz to 200 GPa. The 
recently proposed EOS of Shim et alpll is significantly lower than all the other standards. It is lower than the present analysis by 
10 GPa at V/Vq = 0.7, corresponding to P = 156 GPa. 

Figure shows P(T) along two isochores in comparison with various equations of state. At V /Vq = 0.8, Jamieson et al. 
give somewhat lower P at 300 K with a larger dP/dT, in keeping with their higher y. Heinz and Jeanloz 13 are in the best 
overall agreement with the present pressures, with somewhat lower dP/dT. Shim et alM give a somewhat higher dP/dT with 
P generally low due to their large offset at room temperature. Anderson et alMi adopted the room temperature isotherm from 
Heinz and Jeanloz. 13 Their EOS gives dP/dT substantially too low at compressed volumes. 

We have tabulated P(V, T) in table||l]in the same format adopted by other authors. As discussed, we have extended the range 
of densities to V /Vq = 0.6 and the temperatures to 5000 K, except where such temperatures are thought to lie in the liquid phase. 
Boettger 16 has proposed an extension of the Au 300 K standard to 500 GPa based on LDA calculations. He gives P = 344 GPa 
at V /Vq = 0.6 compared to 329 GPa from our semi-empirical EOS. Given the scatter in the high pressure Hugoniot data, and our 
approximate treatment of melting, this 4.6% difference is probably within the uncertainties of the present analysis at this high 
density. 
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FIG. 6: Room temperature isotherm of Au. Solid curve is present EOS. Solid symbols are data. ■ Open triangles are the EOS of Jamieson 
et alm\ Open squares are EOS of Heinz and Jeanloz 13 Dashed curve is EOS of Shim et aim— 
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FIG. 7: Isochores of Au. Solid curves are the present EOS. Dashed curve is Jamieson et ah- Open squares are Heinz and Jeanloz. Filled 
diamonds are Anderson et aim— Open triangles are Shim et aim—. 
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TABLE II: Tabulated Pressures for Au. Compression T| = 1 — V/V300 where V300 is the volume at T = 300 K and P = 0. Pressures in GPa. 
Values in parenthesis are the first liquid states for each T|, included for interpolation. Remaining liquid states left blank. 

r\ 300 K 500 K 1000 K 1500 K 2000 K 2500 K 3000 K 3500 K 4000 K 4500 K 5000 K 



0.00 


0. 


.00 


1 


44 


5. 


.08 


8 


.73 (12.40) 
















0.02 


3 


.57 


5 


.00 


8. 


.61 


12 


.24 (15.88) 
















0.04 


7, 


.65 


9 


.07 


12. 


.65 


16 


.25 


19.87 (23.50) 














0.06 


12. 


.31 


13 


.71 


17 


.27 


20 


.85 


24.44 (28.04) 














0.08 


17. 


.61 


19 


.01 


22. 


.53 


26 


.09 


29.66 


33.24 (36.83) 












0.10 


23. 


.65 


25 


.03 


28. 


.54 


32 


.07 


35.61 


39.17 


42.74 


(46.33) 










0.12 


30. 


.53 


31 


.90 


35. 


.38 


38 


.88 


42.41 


45.94 


49.50 


(53.06) 










0.14 


38. 


,36 


39 


.71 


43. 


.17 


46 


.65 


50.16 


53.67 


57.20 


60.74 


(64.30) 






0.16 


47. 


.27 


48 


.61 


52. 


.04 


55 


.51 


58.99 


62.48 


65.99 


69.51 


73 


.05 


(76.60) 




0.18 


57. 


.41 


58 


.74 


62 


.15 


65 


.60 


69.06 


72.53 


76.02 


79.53 


83 


.04 


(86.57) 




0.20 


68. 


.96 


70 


.28 


73. 


.67 


77 


.09 


80.54 


83.99 


87.47 


90.95 


94 


.45 


97.96 


(101.49) 


0.22 


82. 


.12 


83 


.43 


86. 


.80 


90 


.21 


93.63 


97.07 


100.53 


104.00 


107 


.48 


110.98 


114.49 


0.24 


97. 


.13 


98 


.43 


101. 


.78 


105 


.17 


108.58 


112.00 


115.44 


118.90 


122 


.37 


125.85 


129.34 


0.26 


114. 


.27 


115 


.55 


118 


.89 


122 


.26 


125.65 


129.07 


132.49 


135.93 


139 


.39 


142.85 


146.34 


0.28 


133. 


.85 


135. 


.12 


138. 


.44 


141 


.80 


145.18 


148.58 


151.99 


155.42 


158 


.86 


162.32 


165.79 


0.30 


156. 


.25 


157. 


.51 


160. 


.81 


164 


.16 


167.53 


170.92 


174.32 


177.74 


181 


.17 


184.62 


188.08 


0.32 


181. 


.92 


183. 


.17 


186. 


.46 


189 


.79 


193.16 


196.53 


199.93 


203.34 


206 


.76 


210.20 


213.65 


0.34 211. 


.38 


212. 


.63 


215. 


.90 


219 


.22 


222.58 


225.95 


229.34 


232.74 


236 


.16 


239.59 


243.04 


0.36 245. 


.26 246. 


.49 


249. 


.75 


253 


.07 


256.42 


259.79 


263.17 


266.57 


269 


.99 


273.41 


276.86 


0.38 


284. 


.30 285. 


.52 


288. 


.76 


292 


.08 


295.43 


298.79 


302.17 


305.57 


308 


.98 


312.41 


315.85 


0.40 


329. 


.37 


330. 


.58 


333. 


.82 


337 


.14 


340.48 


343.85 


347.23 


350.63 


354 


.04 


357.47 


360.91 



V. MELTING AND THE HUGONIOT 



In order to investigate the effects of melting on the Hugoniot, we have constructed a two-phase model free energy. The ion 
free energy in the liquid is based on the assumptions that Cy n = 3Nks, which is reasonable for temperatures near melting, 
and that ASJJ?", the entropy difference between solid and liquid at fixed volume is O.SNIcb- These are empirically based model 
assumptions. 39 A statistical mechanical basis for these observations is discussed by Wallace, 40 who argues for the universality 
of ASy n . Beyond these assumptions, we require the energy of the liquid with respect to the solid to fully determine the liquid 
free energy. We do this by imposing that the melting curve, obtained by equating the pressures and Gibbs free energies of liquid 
and solid, follow the Lindemann rule in the form 

T 

= const, (19) 



coi5(W 2/3 



where V s is the volume on the solidus. The electron excitation free energy is assumed to be the same in the liquid as in the solid. 
Application of this method to Cu leads to a shock melting threshold of 228 GPa, which compares well with 232 GPa obtained 
by Hayes et al. 41 by analyzing sound speed data. The Au Hugoniot has been calculated for the two-phase model allowing for 
coexistence in the shocked state. 42 Figure [8] shows the resulting Hugoniot in the pressure-volume plane. Also shown in the 
figure are the data points and the linear U S (U P ) fit. The dot-dashed curve is the Hugoniot for the solid only. The boundaries 
of the coexistence region are visible as kinks on the solid curve at 280 and 350 GPa. Above 350 GPa, the Hugoniot is in pure 
liquid phase. A similar anomaly is shown for Al in a two-phase calculation by Chisolm et a/4 3 - A significant enhancement of the 
Hugoniot pressure of Au due to melting was also shown in calculations by Godwal et al. 44 , however they predict a much lower 
shock melting threshold than we do. Just above complete melting, the liquid Hugoniot lies above the solid by 30 GPa. The linear 
U S (U P ) curve goes smoothly through the melting region to intersect the liquid data. 

Our two-phase EOS gives good agreement with the high pressure data, which suggests that our cold energy is valid to com- 
pressions of y/Vb « 0.6. Electronic excitations have practically no effect on the P(V) Hugoniot in the solid phase, although 
they significantly affect the temperature. At higher temperatures in the liquid, they act to soften the Hugoniot by absorbing 
energy. Coincidentally, neglecting both electronic excitations and melting gives a Hugoniot that agrees with the linear U S (U P ) 
curve quite well to 650 GPa. This is an accidental cancellation of errors. Ignoring melting and electronic excitations gives a 
temperature that is too high by ~ 10 4 K at this pressure. This offsets the neglect of the pressure enhancement due to melting. It 
is not recommended to ignore either melting or electronic excitations in this high temperature regime. 
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FIG. 8: Hugoniot of Au, including the effects of melting. Solid curve is present two-phase (solid and liquid) EOS. Dot-dashed curve is solid 
EOS only. Dashed curve is linear fit to U s (U p ) data* Solid symbols are measured points^!— 



VI. CONCLUSIONS 



By combining ab initio calculations of elastic moduli and zone boundary phonons with interpolations based on a force constant 
model, we have calculated moments of the phonon frequencies of fee Au to twofold compression. This allows us to calculate the 
associated Griineisen parameters. In particular, we have focused on yo which corresponds to the classical limit. We emphasize 
that the classical limit dominates the thermal pressure in the EOS. We find that the frequently used form y(V) = y(Vo)(V /Vq)^ 
does not represent the volume dependence of y well. We have used an expansion to 2 nd order in V with a physical asymptotic 
limit, which fits the calculated moments well and is consistent with the measured value of y at ambient pressure and temperature. 
Using the theoretical y(V) and electron excitation free energy, we have constructed a semi-empirical free energy for the solid, 
which we believe to be as accurate as possible. This free energy is given as a parameterized closed form expression and the 
resulting P(V, T) is given in tabular form to V/Vb = 0.6. 

Our static lattice energy is empirical, and is verified by comparison with the measured Hugoniot, so that our EOS can be 
regarded as a generalization from the Hugoniot with a physically based y. In the solid phase, electronic excitations have a small 
effect. This, together with the dominance of the classical limit in the vibrational free energy, means that the widely used Mie- 
Griineisen approximation, that (dP/dE) v is independent of T , is accurate. At the highest compressions the physics affecting 
the Hugoniot becomes more complicated. Melting is predicted to begin at 280 GPa and V/Vo = 0.656. Complete melting is 
estimated to lead to a 30 GPa increase in the Hugoniot pressure over the solid. At these densities and higher, the Hugoniot 
temperature is rising rapidly and electronic excitations are playing an increasing role. The highest Hugoniot data are at 580 GPa, 
where the temperature is calculated to be above 2 x 10 4 K. The Mie-Griineisen approximation is not expected to be valid for 
interpolating between the Hugoniot and room temperature at this high compression, because of the strong effects of melting and 
electronic excitations. These effects need to be taken explicitly into account, as has been done here. 

Regarding EOS standards, our analysis gives a room temperature P(V ) curve that agrees well with that of Heinz and Jeanloai^ 
to 200 GPa. Our EOS gives (dP/dT) v generally somewhat larger than theirs, and we believe that in this regard our EOS is 
preferred. The EOS of Jamieson et aim is based on reduction of shock data, and was originally given to 80 GPa. It gives somewhat 
lower pressures than ours, which is due to their use of a Griineisen parameter that is too large at all volumes. Extrapolation of 
their room temperature isotherm to higher pressures 12 is not recommended. Anderson et alM adopted the room temperature 
isotherm of Heinz and Jeanloz 13 , while giving a different thermal pressure. Their EOS gives (dP/dT) v which is substantially 
too low under compression, and is not recommended for high temperatures. Simultaneous compression of Au and Pt showed 12 
that the Pt standard of Holmes et al. 45 gave a pressure higher by 15 and 20 GPa than the Au standards of Heinz and Jeanloz 13 



12 



and Jamieson, respectively, at a pressure of 150 GPa. While the extrapolated Jamieson isotherm is expected to be somewhat 
low, our analysis agrees with the room temperature isotherm of Heinz and Jeanloz, suggesting that the discrepancy between the 
Au and Pt pressures is due to errors in the Pt standard. Given the importance of accurate pressure standards, the status of the Pt 
EOS seems to warrant further investigation. 
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